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Depending on the density reached in the cores of neutron stars, such objects may contain 
stable phases of novel matter found nowhere else in the Universe. This article gives a 
brief overview of these phases of matter and discusses astrophysical constraints on the 
high-density equation of state associated with ultra-dense nuclear matter. 
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1. Introduction 

A forefront area of modern research concerns the exploration of the properties of 
ultra-dense nuclear matter and the determination of the equation of state (EoS)-the 
relation between pressure, temperature and density-of such matter. Experimentally, 
relativistic heavy-ion collision experiments enable physicists to cast a brief glance at 
hot and ultra-dense matter for times as short as about 10 -22 seconds. This is differ- 
ent for neutron stars, which are observed with radio and X-ray telescopes as radio 
pulsars and X-ray pulsars. The matter in the cores of such objects is compressed 
permanently to densities that may be more ten times higher than the densities 
inside atomic nuclei, which make neutron stars natural astrophysical laboratories 
that allow for a wide range of (astro) physical studies and astrophysical phenomena 
(Fig. [U li nked to t he properties of ultra-dense nuclear matter and its associated 
EoS! ^^^^^*^ Of particular interest are neutron stars whose observed prop- 
erties deviate significantly from the norm. Examples of such neutron stars are PSR 
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Fig. 1. The multifaccted connection between high-density nuclear matter and neutron (compact) 
star phenomenal 



J0751+ 1807 whose mass is 2.1 ± 0.2 M©^ neutron star RX J1856. 5-3754 whose 
radius may be > 13 km|H] and XTE J1739-285 whose rotation period may be as 
small as 0.89 msP^l As discussed in this paper, such neutron star data provide an 
excellent opportunity to gain profound insight into the properties of nuclear matter 
at most extreme conditions of density P*^^ 



2. Neutron Star Masses 

Recent results of timing measurements for PSR B1516+02B, located in the globular 
cluster M5, imply a pulsar mass of l-96t^i2 M© (at 68% probability) with a 95% 
probability that the mass of this object is above 1.68 M Q Q2] w hich deviates consid- 
erably from the average masses of binary radio pulsars, Mbrp — 1.35 ±0.04 Mgff^ 
This striking result constrains neutron star masses to at least 1.6 Af Q (2a confi- 
dence level), and even to 1.9 M Q at the la confidence level. The mass and struc- 
ture of spherical, non-rotating neutron stars is calculated by solving the Tolman- 
Oppenheimer-Volkoff (TOV) equation, 

dP(r) _ (s(r) + P(r))[ra(r) + 47rr 3 P(r)) ( ^ 

dr r(r — 2m(r)) 

where m(r) = An J Q r &r' r' 2 e(r') is the gravitational mass inside a sphere of radius 
r. The baryon number enclosed by this sphere is given by 

N(r) = 4tt dr' r' 2 (l - ' n(r') , (2) 

with n(r) the baryon density profile of the star. In order to solve the TOV equa- 
tion one needs to specify the stellar EoS, i.e., the relation between pressure, P, 
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and energy density, e. We apply a b road , modern collection of nuclear equations of 
state, which are compiled in Table lP^Each one f these models was combined at 

Table 1. Collection of nuclear equations of state studied in this paper. The 
entries are: saturation density, n s ; binding energy, ay', incompressibility, K; 
skewness parameter, K'\ symmetry energy, J; symmetry energy derivative, 
L; Dirac effective mass, rarj . 



EoS 


n s 


a v 


K 


K' 


J 


L 


m D 




[fm- 3 ] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


M 


NLp 


0.1459 


-16.062 


203.3 


576.5 


30.8 


83.1 


0.603 


NLp<5 


0.1459 


-16.062 


203.3 


576.5 


31.0 


92.3 


0.603 


D 3 C 


0.1510 


-15.981 


232.5 


-716.8 


31.9 


59.3 


0.541 


DD-F4 


0.1469 


-16.028 


220.4 


1229.2 


32.7 


58.7 


0.556 


KVR 


0.1600 


-15.800 


250.0 


528.8 


28.8 


55.8 


0.805 


KVOR 


0.1600 


-16.000 


275.0 


422.8 


32.9 


73.6 


0.800 


DBHF 


0.1810 


-16.150 


230.0 


507.9 


34.4 


69.4 


0.678 


BBG 


0.1901 


-14.692 


221.6 


-132.4 


36.3 


79.4 




DD-RH 


0.172 


-15.73 


249.0 




34.4 


90.2 


0.686 



sub-nuclear densities with the Baym-Pethick-Sutherland EoSP^The stellar radius, 

R, is defined as that stellar location where pressure vanishes, P(R) = 0. The star's 

gravitational mass and its total baryon number are thus given by M = m(R) and 

N = N(R), respectively. Most of the models sh own in Tab le 1 are derived in the 

framework of the relativistic mean-fiel d ap proach! 1^9 | 20 | 21 | a Uowing for non- linear 

fool 

(NL) self-interactions of the a meson For model NLp, the isovector part of the 
interaction is described entirely in terms of p meson exchange. This is different for 
NLpi5 where the isovector part of the interaction is described in terms of both p and 
5-meson exchange. The latter is generally neglected in RMF modelsPSl RMF mod- 
els with density dependent input parameters (coupling constants and masses) are 
represented in Table 1 by four different models from two classes, where in the first 
one density dependent meson couplings are modeled such that several properties of 
finite nuclei (binding energies, charge and diffraction radii, surface thicknesses, neu- 
tron skin in 208 Pb, spin-orbit splittings) can be fitted!^ D 3 C additionally contains 
a derivative coupling which leads to momentum-dependent nucleon self-energies, 
and DD-F4 is modeled such that the flow constraint from heavy-ion collisions is 
fulfilledPS] The second class of these models is motivated by the Brown- Rho scaling 
assumptiorJ2Sl that not only the nucleon mass but also the meson masses should 

|97| 

decrease with increasing density. In the KVR and KVOR models^- 1 these depen- 
dences are related to a nonlinear scaling function of the c-meson field such that the 
EoS of symmetric nuclear matter and pure neutron matter below four times the 
saturation density coincide with those of the Urbana-Argonne groupPSlln this way 
the latter approach builds a bridge between the phenomenological RMF models and 
a microscopic EoS built on realistic nucleon-nucleon forces. The RMF models are 
contrasted with several variational models for the EoS such as APR^, WFF^, 
ppcpS a relativistic Dirac-Brueckner-Hartree-Fock (DBHF) modef^, and a non- 
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Fig. 2. Mass versus central density of neutron stars, computed from Eq. JTJ for EoS shown in 
Table 1. Crosses refer to the maximum-mass star of each sequence, filled dots mark the critical 
masses and central densities beyond which the direct Urea (DU) cooling process becomes possible. 

relativistic Brueckner-Bethe-Goldstone (BBG) modelP^ As shown in Fig. [2j none 
of these values falls below the 2a mass limit of 1.68 M Q for PSR B1516+02B, and 
even at the la mass limit of 1.84 M@ the softest EoS NLp and NLpS cannot be 
excluded. We point out that if a pulsar with a mass exceeding 1.8 — 1.9 Mq at the 
2a or even 3a level would be observed in the future, this would impose severe con- 
straints on the stiffness of the nuclear EoS. For the set of EoS tested here, only the 
stiffest models, i.e. D 3 C, DD-F4, BBG, and DBHF would remain viable candidates. 

3. Gravitational Mass Baryon Number Relation 

It has been suggested that pulsar B in the double pulsar system J0737-3039 may 

IOOI 

serve to test models proposed for the EoS of superdense nuclear matter!^ One of 
the interesting cha racteri stics of this system is that the mass of pulsar B is merely 
1.249 ± 0.001 M Q .12^1351 Such a low 

mass could be an indication that pulsar B 
did not form in a type-II supernova, triggered by a collapsing iron core, but in 
a type-I supernova of an ONeMg white dwarf driven hydrostatically unstable by 
electron captures onto Mg and NeP^ The well-established critical density at which 
the collapse of such stars sets in is 4.5 x 10 9 g/cm 3 . Assuming that the loss of matter 
during the formation of the neutron star is negligible, a predicted baryon mass for 
the neutron star of Mat = 1.366 — 1.375 Mq was derived in Ref.1331 This theoretically 
inferred baryon number range together with the star's observed gravitational mass 
of M = 1.249 ± 0.001 Mq may represent a most valuable constraint on the EoS, 
provid ed the above key assumption for the formation mechanism of the pulsar B is 
correct 133136] If 

so, any viable EoS proposed for neutron star matter must predict 
a baryon number in the range 1.366 ^ Af/v ~ 1.375 Mq for a neutron star whose 
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Fig. 3. Gravitational mass— baryon mass relationship!^! for PSR J0737-3039 ( B) co mpared to a 
set of rclativistic EoS and four nonrelativistic EoS (BBG, APR, FPS and WFF)!Sl 

gravitational mass is M — 1.249 ± 0.001 Mq. We may contrast this result with 
that of an independent calculation}^ where for pulsar B a baryon mass of = 
1.360 ±0.002 Mq has been obtained. The authors of Ref . 1531 discussed caveats such 
as baryon loss and variations of the critical mass due to carbon flashes during the 
collapse. The effect of 1% and 2% mass loss on the usefulness of this constraint to 
exclude model EoS is shown in Fig. [3] 

4. Mass-Radius Constraints from Neutron Stars in LMXBs 

Aside from neutron star masses, kilohcrtz quasi-periodic brightness oscillations 
(QPOs) seen from more than 25 neutron star X-ray binaries (LXMBs) can be used 
to put additional constraints on the high-density EoS. A pair of such QPOs is often 
seen from these systems!^ In all currently viable models for these QPOs, the higher 
QPO frequency is close to the orbital frequency at some special radius. For such 
a QPO to last the required many cycles (up to ~ 100 in some sources), the orbit 
must be outside the star. According to general relativity theory the orbit must also 
be outside the innermost stable circular orbit (ISCO). Gas or particles inside the 
ISCO woul d spira l rapidly into the star, preventing the production of sharp QPOs. 
This impliefP^QI 

that the observation of a source whose maximum QPO frequency 
is t'max limits the stellar mass and radius to 

M < 2.2 M Q w °°^ z (1 + 0.75j) , R < 19.5 km 10 ° m ° m Hz (l + 0.2 j) . (3) 

The quantity j = cJ/GM 2 (with J the stellar angular momentum) is the dimcn- 
sionless spin parameter, which is typically in the range between 0.10 and 0.2 for 
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Fig. 4. Left panel: Mass-radius constraints from thermal radiation of the isolated neutron star 
RX J1856 (grey hatched region) and from QPOs in the LMXBs 4U 0614+09 (wedge-like green- 
hatched area) and 4U 1636-536 (orange hatched region). For 4U 1636-536 a mass of 2.0 ±0.1 Mq 
is obtained so that the weak QPO constraint would exclude the NLp and NLp<5 EoS, whereas 
the strong one renders only the stiffest EoS D 3 C, DD-F4, BBG and DBHF viable. XTE J1739- 
285 would favor a soft EoS, which is controversial, however (see text). Right panel: Same as left 
panel, but for hybrid stars with a DBHF hadronic shell and a 2SC color superconducting NJL 
quark matter core computed for different coupling strengths, f)jy and rjy. The dotted lines display 
surface redshifts, z, ranging from 0.1, ... , 0.7. 



these systems. Equation © implies that for given observed value i/ max the mass 
and radius of that source must be inside a wedge-shaped shown in Fig. [5] 

Since wedge becomes smaller for higher z/ max , the highest frequency ever observed, 
1330 Hz for 4U 0614+09lS2, places the strongest constraint on the EoS. As can be 
seen from Fig. [U the current QPO constraints do not rule out any of the EoS con- 
sidered here. However, because higher frequencies imply smaller wedges, the future 
observation of a QPO with a frequency in the range of ~ 1500 — 1600 Hz would rule 
out the stiffest of our EoS. 

If there is evidence for a particular source that a given frequency is close to the 
orbital frequency at the ISCO, then the mass is known to a good accuracy, with 
uncertainties arising from the spin parameter. This was first claimed for 4U 1820- 
30 but complexities in the source phenomenology have made this controversial. 
More recently, careful analysis of Rossi X-ray Timing Explorer data for 4U 1636- 
536 and other sourced^ has suggested that sharp and reproducible changes in 
QPO properties are related to the ISCO. If so, this implies that several neutron 
stars in low-mass X- ray binaries should have gravitational masses between 1.9 Mq 
and possibly 2.1 M F^I i n Fig. [4] we show the estimated mass of 2.0 ± 0.1 M© for 
4U 1636-536. 

Recently, mass-radius constraints have been reported for the accreting compact 
stars XTE J1739-285 44 and SAX J1808.4-365cPI (SAX J1808 for short) which are 
based on the identification of the burst oscillation frequency with the spin frequency 
of the compact star. It is almost impossible to fulfill the constraints from RX J1856 
(and other high-mass candidates) and SAX J1808 (which favors a soft EoS) simul- 
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taneously so that the status of SAX J1808 is currently controversial. It is likely that 
the small radius estimate of Leahy et al@3 is a consequence of the underestimation 
of higher harmonics when only timing data are analyzed and not also the energy 
spectra!^ Neutron star XTE J 1739-285 and its possible implications for the EoS 
are discussed in more detail below. 

5. Mass-Radius Relation Constraint from RX J1856. 5-3754 

The nearby isolated neutron star RX J1856. 5-3754 (RX J1856 for short) belongs to 
a group of seven objects which show a purely thermal spectrum in X-rays and in 
optical- UV. This allows the determination of R^/d, the ratio of the photospheric 
radius to the distance d of the object, if the radiative properties of its photo- 
sphere are known. RX J1856 is the only object of this group which has a measured 
distance obtained by Hubble Space Telescope (HST) astrometry. After the distance 
of 117 pc^^ became known several groups pointed out that the blackbody radius 
of this star is as large as 15 to 17 km. Although both the X-ray and the optical-UV 
spectra are extremely well represented by blackbody functions they require different 
emission areas, a smaller hot spot and a larger cooler region. The overall spectrum 
could also be fitted by blackbody emission from a surface showing a continuous 
temperature distribution between a hot pole and a cool equator, as expected for a 
magnetized neutron star. The resulting blackbody radii are 17 km (two blackbodies) 
and 16.8 km (continuous temperature distribution jH| i n t n j s paper we adopt the 
result of the continuous temperature fit, Roo = 16.8 km. More recent HST observa- 
tions of RX J1856 indicated larger distances of up to 178 pc. A distance of around 
140 pc for RX J1856 is considered a conservative lower limit.HHI For a distance of 
140 pc the corresponding radius is 17 km!^ Although some questions-in particular 
that of the distance-are not yet finally settled, the recent data support an unusually 
large radius for RX J1856. 5-3754. 

6. Surface Redshift 

An additional test to the mass-radius relationship is provided by measurements of 
the gravitational redshift of line emissions from the surfaces of neutron stars. The 
gravitational redshift, z is given in terms of the star's gravitational mass and radius 
according to 



In the right panel of Fig. [5l the dependence of z on star mass is shown fo r two 
representative EoS. The disputed measurement of z = 0.35 for EXO 0748-676^112 
allows for both an unconfined hadronic core composition as well as a deconfined 
quark-core composition.-^ A measurement of z > 0.5 could not be reconciled with 
the hybrid star models suggested here, while the hadronic model would not be 
invalidated by redshift measurements up to z = 0.6. 




(4) 
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Fig. 5. Left panel: Moment of inertia versus mass. The highlighted mass regions correspond to 
the double pulsar J0737+3039 A+B and the pulsar B1516+02BA. Right panel: Surface redshift 
versus gravitational mass for a purely hadronic star (DBHF) as well as hybrid stars with 2SC 
color superconducting NJL q uark m atter cores. The dotted line labeled EXO 0748-676 displays 
the disputed redshift z = 0.35.l 50 l 51 l 



7. Moment of Inertia 

The observed data of the relativistic double puls ar PSR J0737+3039 did allow a 
determination of the moment of inertia (Molji^J f th at star , which may be used 
to put further constraints on the EoS of neutron stars Results for the Mol 
computed for the EoS of this paper are shown in Fig. [5] Due to the fact that the 
mass 1.338 Mq of PSR J0737+3039 A is in the vicinity of the suggested critical 
mass region, the quark matter core is small and the expected Mol of the hybrid star 
will be practically indistinguishable from that of a pure hadronic one. The situation 
would improve if the Mol could be measured for more massive objects, because the 
difference in the Mol of both alternative models for masses as high as 2 Mq could 
reach the 10% accuracy level. 



8. Rotational Frequencies 

An absolute upper limit on the spin frequency of a pulsar is given by the mass 
shedding limit, at which the velocity of the stellar surface equals that of an orbit- 
ing particle suspended just above the surface. For a rigid Newtonian sphere this 
frequency is given by the Keplerian frequency)^ 

v K = (27r)- 1 V GM/R 3 = 1833 (M/Af Q ) 1/2 (i?/10 km)" 3/2 Hz . (5) 

This formula was found to describe the mass shedding points for a sample of neu- 
tron star EoS extremely well However, both deformation and general relativistic 
effects are very important so that Eq. ([5]) needs to be modified. It has been founcP^ 
that with a coefficient of 1045 Hz, Eq. (JSJ) approximately describes the maximum 
rotation rate for a star of mass M and non-rotating radius R, independently of 
the EoS. The observation of rapidly rotating pulsars can therefore constrain the 
compactness and might eventually lead to the elimination of EoS that are too stiff, 
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since the latter lead to stars that are too big. If the recent discovery of burst os- 
cillations with a frequency of 1122 Hz in the X-ray binary XTE J1739-28512I and 
their identification with the star's spin frequency turns out to be correct, this object 
would spin at a rotational period of 0.89 ms, rendering it the first sub-millisecond 
pulsar ever observed, imposing severe constraints on the EoS. The boundary for the 
mass-radius relationship reads in this case^l 

R < 9.52 (M/M Q ) 1/3 km . (6) 

As can be seen from the left panel in Fig. [H this constraint would rule out all neu- 
tron stars of masses M < 1.75 M Q . Only neutron stars computed for the stiffest 
EoS of our collection (D 3 C, DD-F4, BBG and DBHF) remain viable. In conclud- 
ing this section, we mention that an important test of the transport prop erti es of 
the matter inside sub- millisecond pulsars arises from the r-mode instability!^ This 
constraint on superconducting quark matter in neutron stars has recently been 
discussed for the putative sub-millisecond pulsar XTE J1739-285pHl where it was 
argued that this object must be either a strange star or a quark-hadron hybrid 
star. In passing we mention that rotational instabilities in rotating stars, known 
as gravitational radiation driven instabilities, set a m ore stringent limit on rapid 
stellar rotation than mass shedding from the equatorPUl These instabilities orig- 
inate from counter-rotating surface vibrational modes which at sufficiently high 
rotational star frequencies are dragged forward. In this case gravitational radiation, 
which inevitably accompanies the aspherical transport of matter, does not damp 
the instability modes but rather drives them. Viscosity plays the important role of 
damping these instabilities at a sufficiently reduced rotational frequency such that 
the viscous damping rate and power in gravity waves are comparable. The most 
critical instability modes that are driven unstable by gravitational radiation are the 
/-modes and the recently discovered r-modes. 1 ^ The latter may severely const rain 
the composition of compact stars that would rotate at sub-millisecond periodsPSl 



9. Mass Clustering 

Compact stars in binary systems in general undergo during their evolution a stage 
with disc accretion leading to both, spin-up and mass increase. Initial indications 
for a spin frequency clustering in low-mass X-ray binary systems, reported by mea- 
surements with the Rossi-XTE, have lead to the suggestion to interprete such a cor- 
relation as a waiting-point phenomenon where star configurations cross the border 
between pure neutron stars and hybrid stars in the spin frequency-mass planeP^ 
A systematic analysis of the critical line for a deconfinement phase transition in 
the phase diagram for accreting compact starJS^I has revealed that the suggested 
population clustering due to the phase transition shall rather lead to a m ass cluster- 
ing effect. For strange stars, however, such an effect shall be absent. For generic 
polytropic forms of the EoS of quark and hadronic matter, the relationship between 
softness or hardness of the EoS and the structure of this phase diagram has been 
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Fig. 6. Rotational frequency, Q, versus mass, M, of neutron stars. Stable stars are located in 
the region bordered by the maximum rotation frequency f! max (M), where mass-shedding from 
the equator occurs, and the maximum mass M max (f2) at which gravitational instability against 
collapse to a black hole sets in. The dash-dotted (dash-double-dotted) line marks the onset of the 
generation of a quark matter core in the star's interior for a diquark coupling parameter r)r> = 1.020 
(nrj = 1.017), see Fig. [4] for the corresponding mass-radius relations of the non-rotating stellar 
sequences. The data with error bars correspond to binary radio pulsars (black dots) and neutron 
stars with a white dwarf binary (green squares). The correlation of the distribution of objects with 
the line for a deconfinement phase transition could be interpreted as a wa iting point phenomenon 
during the accretion evolution of the compact stars (mass clustering) .ESI 



demonstrated in Ref. [63] For the specific example of the DBHF hadronic EoS and a 
color superconducting stiff quark matter EoS, we show the phase diagram of rotating 
compact stars in Fig. [6l For curiosity we show also the masses and spin frequencis 
for compact stars in binary radio pulsars and neutron star-white dwarf binaries 
and observe an interesting correlation of the distribution of these objects with the 
critical deconfinement phase transition line for diquark coupling r/D = 1.02. There 
might be other reasons for the mass clustering of spinning pulsarsp^ but the sug- 
gestion to relate it to a phase transition in the interior can not be excluded. If true, 
it would be a strong constraint for the hybrid EoS and suggest a critical density for 
deconfinement in compact stars at about 0.4 fm -3 (i.e., 2.5 n s ). We want to point 
out that the non-zero quark chemical potential (fi ^ 0) domain is a rather poorly 
understood region of the QCD phase diagram. Fundamental approaches, lik e solving 
the in-medium QCD Schwinger-Dyson equations in a specific QCD model E 5 | 69 | 70 | 
to obtain a quark matter EoS are demanding. This well justifies the description of 
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Fig. 7. Hadronic star cooling curves for DBHF (left) and DDF-4 (right) model EoS. Different 
lines correspond to compact star mass values indicated in the legend (in units of Mq), data points 
with error bars are taken from Ref. EH. The hatched trapeze-like region represents the brightness 
constraint (BC)!S5]The onset of the DU process entails a high sensitivity of TA curves to small 
mass variations known as DU problerrJSSl 

a quark matter phase within three-flavor NJL- type model s, giving access even to 
diquark pairing channels at reasonable expense J2XE2E3I73I 



10. Stellar Cooling 

Neutron stars cool for the first 10 5 years through neutrino emission from their cores, 
which is followed at later times by photon emission. We simulate this behavior by 
using the cooling code developed in Refs. [75l [TBI The cooling calculations are per- 
formed for different hadronic EoS, different assumptions about superconductivity, 
and different stellar crust models. Moreover, we try to use consistent microscopic 
inputs. (For systematic field-theoretical approaches to the neutrino cooling problem 
in neutron stars see, for instance, Ref. [77J [5l [78j US) The main neutrino cooling 
processes in hadronic matter are the direct Urea (DU), the medium modified Urea 
(MMU) and the pair breaking and formation (PBF). For quark matter, the main 
cooling processes are the quark direct Urea (QDU), quark modified Urea (QMU), 
quark bremsstrahlung (QB) and quark pair formation and breaking (QPFBpQl 
Also the electron bremsstrahlung (EB) , and the massive gluon-photon decay^ are 
included. 

The 1 Sq neutron and proton gaps in the hadronic shell are taken according to the 
calculations in Ref. [82] corresponding to the thick lines in Fig. 5 of Ref. l75l However, 
the 3 P2 gap is suppressed by a factor 10 compared to the BCS model calculation 
in Ref. 82, consistent with arguments from a renormalization group treatment of 
nuclear pairing!^ Without such a suppression of the z Pi gap the hadronic cooling 
scenario would not fulfill the TA constraint!^ 

The possibilities of pion condensation and of other so called exotic processes are 
included in the calculations for purely hadronic stars but do not occur in the hybrid 
star calculations since the critical density for pion condensation exceeds that for 



August 10, 2008 17:43 WSPC/INSTRUCTION FILE paper 



12 D. Blaschke, T. Klaehn, and F. Weber 




Fig. 8. Density dependence of the pairing gaps in nuclear matter together with that of the 
hypothetical X-gap in quark matter (left). Mass-central density relation for the two hadronic EoS 
models DBHF and DD-F4 (right). The dot indicates the onset of the DU process. 



deconfinement in our caseM-^ While the hadronic DU process occurs in the DBHF 
model EoS for all neutron stars with masses above 1.27 Mq, it is not present at 
all in the DD-F4 model, see the right panel of Fig. [H We account for the specific 
heat and the heat conductivity of all particle species whose presence is predicted by 
/3-equilibrium. In addition, in the case of quark matter the contributions of massless 
and massive gluon-photon modes are taken into account too. 

In the 2SC phase only the contributions of quarks forming Cooper pairs (say 
red and green) are suppressed via huge diquark gaps, while those of the remaining 
unpaired blue color lead to a so fast cooling that the hybrid cooling scenario becomes 
unfavorable!^ We thus assume the existence of a weak pairing channel such that 
in the dispersion relation of hitherto unpaired blue quarks a small residual gap 
can appear. We call this gap Ax and show that for a successful description of the 
cooling scenario Ax has to have a density dependence. We have studied the ansatz 
Ax = Ao exp [— a(/i//i c — 1)], where /i is the quark chemical potential, /i c = 330 
MeV. For the analyses of possible models we vary the values of a and Ao, given in 
the Table 1 of Ref. 86] and shown in the left panel of Fig. [8] 



11. Temperature-Age (TA) Test 

We consider the cooling evolution of young neutron stars with ages t ~ 10 3 — 10 6 yr 
which is governed by the emission of neutrinos from the interior for t 5; 10 5 yr 
and thermal photon emission for t > 10 5 yr. The internal temperature is on the 
order of T ~ 1 keV. This is much smaller than the neutrino opacity temperature 
T op ac ~ 1 MeV as well as the critical temperatures for superconductivity in nuclear 
(T c ~ 1 MeV) or quark matter (T c ~ 1 — 100 MeV). Therefore, the neutrinos are 
not trapped and the matter is in a superconducting sta te. In Fig. [8] we show the 
density dependence of the pairing gaps in nuclear m atteJ^^ together with those 
of the hypothetical X-gap in quark matter JHZE5IHS1 The phase transition occurs at 



August 10, 2008 17:43 WSPC/INSTRUCTION FILE paper 



Constraints on the nuclear EoS from neutron star observables 13 




Fig. 9. Cooling curves for hybrid star configurations with 2SC+X pairing pattern and X-gap 
model I (left) versus model IV (right). For the gaps see the left panel of Fig. [7] The grey value for 
the shading of the mass bin areas corresponds to the probability for that mass bin value in the 
population synthesis model of Ref. 1881 for details see Ref. 1661 



40- 
20- 

40- 
20- 



40- 
20-, 



DBHF 



DDF-4 



Popov et al. (2006) 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 
M [M. J 



Popov et al. (2006) 



1.2 1.4 1.6 
M [M ] 



Fig. 10. NS mass spectra extracted from the distribution of cooling data for both hadronic EoS 
models (left panel) and for hybrid stars with X-gap models I- IV (right panel). For co mpar ison, the 
mass distribution of young, nearby NS from the population synthesis of Popov et al. 12^1 i s shown 
at the bottom of the panels. 



the critical density n c = 2.75 uq = 0.44 fm -3 . 

In Fig. [7] we present temperature-age (TA) diagrams for two different hadronic 
models. Figure [5] shows the TA diagram for two hybrid star cooling models which 
are presented in Ref. [86l We note that the T m — T s relationship between the tem- 
peratures of the inner crust and th e st ellar surface has been chosen according to 
a Tsuruta's formula, for details see ^1. The "TA test" is fulfilled when each data 
point ought to be explained with a cooling curve of an admissible configuration. The 
TA data points are taken from Ref. [64] The hatched trapeze-like region represents 
the brightness constraint (Bcj^. For each model nine cooling curves are shown 
for configurations with mass values corresponding to the binning of the population 
synthesis calculations explained in Ref. [86l A so called logN-LogS distribution con- 
straint has been considered in Ref. [86] for the hybrid cooling scenario above, and it 
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has been suggested to use the marking of TA diagram with five grey values in order 
to encode the likelihood that stars in that mass interval can be found in the solar 
neighborhood, in accordance with the population synthesis scenario, see Fig. [9] The 
darkest grey value, for example, corresponds to the most populated mass interval 

I. 35 to 1.45 Mq predicted by the mass spectrum used in population synthesis. It has 
been suggested in Refl66l that a mass spectrum can be derived from the sequences 
of TA curves of a given compact star cooling theory, see Fig. [TUJ This example 
demonstrates how in principle cooling simulations together with data from obser- 
vation and population synthesis simul atio ns could discriminate pairing patterns for 
quark matter phases. A recent review^2l gives a flavor of the fascinating topics in 
the discussion of color superconducting phases in compact stars. Many new devel- 
opments in the theory and observations of compact stars like, e.g., the discussion 
of single flavoJ^ and strange quark matteJ^ phases for providing possible deep 
crustal heating mechanis ms to e xplain the puzzling phenomena of superbursts and 
cooling of X-ray transient d^^ will deepen our understanding of the high-density 
nuclear EoS. 
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